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NON-LINEAR  TRANSVERSE  ELECTRON  BEAM  DYNAMICS 
IN  A  MODIFIED  BETATRON  ACCELERATOR 


I.  Introduce ion 

High  energy  accelerators  capable  of  producing  high  current  electron  beams 
are  rapidly  becoming  an  active  area  of  research.  The  motivation  for 
developing  these  devices  is  related  to  potential  applications  of  high  current 

beams  to  the  generation  of  high  power  coherent  radiation1,  x-ray  radiography 

2 

and  national  defense  « 

Among  the  various  accelerating  schemes  that  have  the  potential  to  produce 
ultra-high  power  electron  beams,  Induction  accelerators^  appear  to  be  the  most 
promising.  Induction  accelerators  are  inherently  low  Impedance  devices  and 

V 

thus  are  Ideally  suited  to  drive  high  current  beams.  The  acceleration  process 
is  based  on  the  inductive  electric  field  produced  by  a  time  varying  magnetic 
field.  The  electric  field  can  be  either  continuous  or  localized  along  the 
acceleration  path. 

Quite  naturally,  induction  accelerators  are  divided  into  linear  and 

cyclic.  The  linear  devices  are  in  turn  divided  into  Astron-type*-8,  Radlac- 
9  10  11  12 

type  *  and  auto-accelerator  *  .  In  the  first  type,  ferromagnetic 

induction  cores  are  used  to  generate  the  accelerating  field,  while  "air  core” 

cavities  are  used  in  the  second.  In  the  auto-accelerator  the  air  core 

cavities  are  excited  by  the  beam's  self  fields  rather  than  external  fields. 

13—15 

Similarly,  cyclic  devices  can  be  divided  into  conventional  and  modified 

betatrons18”^®.  The  field  configuration  in  the  modified  betatron  Includes,  in 

addition  to  the  time  varying  betatron  magnetic  field,  which  is  responsible  for 

the  acceleration,  a  strong  toroidal  magnetic  field  that  substantially  Improves 

the  stability  of  the  accelerated  beam. 

The  linear  dynamics  of  high  current  electron  rings  in  modified  betatron 

fields,  with  and  without  stellarator  fields,  has  been  studied  extensively 
Mnuacript  approved  July  26, 1984. 
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during  Che  last  few  years  .  These  studies  are  based  on  the  linearized 
equations  of  motion,  i.e.,  they  assume  that  the  electron  ring  is  confined  near 
the  minor  axis  of  the  torus. 

In  this  paper  we  study  the  transverse  ring  dynamics  in  a  modified 


betatron  accelerator  using  a  different  approach.  The  ring  orbits  are  not 
determined  from  the  equations  of  motion  but  rather  from  the  two  exact 
constants  of  the  motion  and  the  potentials  at  the  centroid  of  the  ring.  The 
main  advantage  of  the  present  approach  is  that  the  ring  orbits  can  be 
determined  over  the  entire  minor  cross-section  of  the  torus  and  not  only  near 
its  minor  axis.  The  topology  of  orbits  near  the  wall  of  the  toroidal  vacuum 
chamber  is  of  vital  Importance  during  Injection,  since  optimization  of  the 
confining  region  requires  the  beam  to  be  Injected  far  away  from  the  minor  axis 
of  the  torus^’^.  It  has  been  found  that  the  shape  of  the  ring  orbits,  in 
the  transverse  to  the  toroidal  magnetic  field  plane,  can  be  very  complex,  in 
particular  in  the  high  current  limit. 

II.  Constants  of  the  Motion 

Consider  an  electron  ring  inside  a  perfectly  conducting  torus  of  circular 

cross  section  as  shown  in  Fig.  1.  The  center  of  the  ring  is  located  at  a 

2 

distance  Ar.Az  from  the  minor  axis  of  the  torus.  The  kinetic  energy  me  y  of  a 
reference  electron  that  is  located  at  the  position  r,z  varies  according  to  the 
equation 


i 

V.V 


nc2  (r,z)  -  -  |m|  v  •  I  (r,z),  (1) 

where  E  (r,z)  is  the  total  electric  field  at  the  position  of  the  reference 
electron.  The  electric  field  is  related  to  the  space  charge  t  and  magnetic 
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vector  potential  A  by 


E  (r,z) 


_  7«  _  ± 

c  3t  * 


where  the  total  time  derivative  of  ♦  is  given  by 


d$  3$  ♦  _ 

dj  -  -jj  ♦  V  •  7». 


For  the  problem  of  interest,  the  accelerating  fields  vary  slowly  in  time  and 


ii.il  .  0 

3t  3t 


Combining  Eqs.  (1)  to  (4),  we  obtain 


dv(r.z) 


-H 


y(r,z)  -  -I— jr  $(r,z)  -  constant. 


According  to  Eq.  (5)  the  sum  of  the  kinetic  and  potential  energy  of  the 


reference  electron  is  conserved. 


Since  the  fields  of  the  modified  betatron  configuration  are  Independent 


of  the  toroidal  angle  8,  the  canonical  angular  momentum  Pg  is  also  a  constant 


of  the  motion,  i.e., 


Pg  -  y  mr  Vg  -  -i|i  r  Ag  ■  constant, 


where  Ag  is  Che  toroidal  component  of  the  total  magnetic  vector  potential 
and  vQ  is  the  toroidal  velocity  of  the  reference  electron. 

Assuming  that  »v  and  eliminating  y  from  Eqs.  (5)  and  (6),  it  is  obtained 

[(  Ag(r,z) )  +l]  - -U|-  *(r,z)  -  constant,  (7a) 


or,  at  the  centroid  of  the  ring 


P.  |  |  2  V2  || 

^mcR  +  A0<R»Z))  +1]  "  •CR.Z)  -  constant. 


For  very  high  energy  beams,  l.e.,  when  y  »  1,  Eq.  (7b)  is  reduced  to 


—g  +  [Aq(R,Z)  -  «(R,Z)]  -  constant. 


This  non-linear  conservation  law  can  furnish  very  useful  information  on 
the  motion  of  the  ring  in  the  r,z  plane,  provided  that  the  potentials  Ag 
and  9  at  the  center  of  the  ring  are  known.  It  should  be  noticed  that  Eqs.  (7) 
are  Independent  of  the  toroidal  magnetic  field. 


III.  The  Potentials 


In  Eq.  (7),  the  total  magnetic  vector  potential  Ag(r,z)  is 


Afl(r,z)  -  Ao*t(r,z)  +  Aflself(r,z), 


where  AQ*xt(r,z)  is  the  external  and  AgSe^(r,z)  is  the  self  magnetic  vector 
potential . 


It  is  assumed  that  Che  becatron  magnetic  field  is  described  by 


,  ext,  N 
Ag  ( r  j  z) 


B 

zo 


(8) 


where  Bzo  is  the  magnetic  field  at  r-rQ,  z  -  o  and  n  is  the  external  field 
index,  i.e.. 


n 


zo 


3B, 

(tt). 
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For  a  cylindrical  electron  beam  Inside  a  straight  perfectly  conducting 
cylindrical  pipe  of  circular  cross  section,  the  self  potentials  can  be 
computed  exactly,  even  for  large  beam  displacements  from  the  minor  axis  of  the 
torus.  In  the  local  coordinate  system  p,$  the  self  potentials  inside  the 
beam,  i.e.,  for  |p-£|  <  rb  are  given  by 

A8eelf<P,4.)  -  -  2 1 e |  ■  09  (V  in  ~  -  tpV-2p  ,4  c™  {4~«)1 

rb  2rb 

“  2  if)  (“)  l~l  cos  (♦-a)},  (9a) 

*  a 


and 


♦(.,♦)  -  -  2 |e |  »,  {  V2  *  tn  ^  '  2  S.  *  <»--)) 

“  2  (f)  (“)  a"1  cos  (♦><*)},  (9b) 
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At  the  beam  center,  l.e.,  for  p  *  A  and  $  *  a,  Eqs.  (9a)  and  (9b)  become 


Agelf(R,Z)  -  -2  |e|N£  g0  {  V2  +  An  ^  +  An  [l  - 


(R-r  )2  +  Z2 


]},  (10a) 


(R-r)2  +  Z2 

*(R,Z)  -  -  2  | e |  {V2+  In  *  +  An  [l - 2__ -  ]}, 

b  a 


(10b) 


where  is  the  linear  electron  density,  r^  is  the  minor  radius  of  the  beam, 
a  is  the  minor  radius  of  the  conducting  pipe  and  00  -  vd/c* 

Figure  2  shows  the  ratio  -A0e^/2N£  |e|  (v0^)  from  Eq.  (10a)  at  Z  -  o, 
together  with  results  from  the  computer  code  PANDIRA.  This  code  solves  the 
differential  equations  for  X  and  *  in  a  non-uniform  triangular  mesh  in  r-z 
coordinates  and  its  present  version  has  been  developed  by  R.F.  Holsinger.  The 
various  parameters  for  the  runs  shown  in  Figs.  2  to  5  are  listed  in  Table  1. 
The  agreement  between  Eq.  (10a)  and  the  numerical  results  is  excellent.  The 
maximum  difference  between  the  analytical  and  numerical  results  is  less  than 
0.4%. 

self  self 

Figure  3  shows  the  stream  function  ♦  -  RA0  .  In  contrast  to  A0  ,  the 
stream  function  i|»  peaks  away  from  the  minor  axis  of  the  torus.  The  radial 
displacement  of  the  peak  can  be  computed  from  8i|>/3R  *  0  and  is  given  by  the 


relation 


A  -  (a/2)  (a/r  )  (V2+  An  a/r  ). 


This  relation  predicts  that  i>  peaks  2.8  cm  from  the  minor  axis,  which  is  not 
by  far-off  from  the  2.6  cm  of  Fig.  3. 

Results  for  the  electrostatic  potential  are  given  in  Fig.  4.  Again  the 
agreement  between  Eq.  (10b)  and  numerical  results  is  excellent.  The  maximum 


I 


difference  A$  between  the  electrostatic  potential  computed  from  PANDIRA  and 
that  of  Eq.  (10b)  is  less  than  2Z. 

Figure  5  shows  the  stream  function  <(*  for  a  torus  with  a  major  radius  rQ  * 
32  cm.  In  agreement  with  the  approximate  expression  for  Ap,  the  displacement 
of  the  peak  increased  by  about  a  factor  of  3. 

To  obtain  a  better  understanding  of  the  potentials  inside  a  perfectly 
conducting  torus,  we  solved2^  the  differential  equations  for  $  and  X  to  first 
order  in  the  ratio  a/R,  but  to  any  order  in  the  normalized  displacement  A/a. 

For  a  constant  particle  density  nQ  ring  and  to  second  order  in  A/a,  the 
electrostatic  potential  and  the  stream  function  \Ji  at  the  center  of  the  ring 
are  given  by 


(R-r  )2  +  Z2 

4(R,Z)  -  -  2  N£  | e |  [V2+  in  (a/rb) - ^ - 

& 


(He) 


(R-r  )2  +  Z2 

<KR,Z)  ■  -  2  N4  |e|  R  8g[  V2  +  in  (a/rb) - - - 

aT 

3r2  (-R-r  > 

- r-  ^”5 — Jj,  for  0  -  constant  (lib) 

8a2  R 


(R— r  )2  +  Z2  r  2  (R-r  ) 

iKR.Z)  »  -  2N4|e|R08  [1/2  +  in(a/rb) - 2- - - g2-  j, 

a  8a 


for  Jq  ■  constant 


(11c) 


Similarly,  the  fields  at  the  centroid  of  the  ring  are  given  by 


2  eN«  j  V, 

a  >•  a  *n  rb  +  8Ra  ^ 


2  |  e  |  N 


B 


2HV9  2 

a  ''a-'* 


2  e  N,Ss  (R-r  ) 

B - 4-2  I - 2_ 

z  a  L  a 


-(#  t1”  -r^) 

D 


+  (3rb2/8Ra)J,  for  9  -  constant. 


(12a) 


(12b) 


(12c) 


( 12d) 


and 


( 12e) 


The  toroidal  term  in  Eq.  (11)  is  very  small  for  the  parameters  of  Interest 
and  therefore  it  is  not  surprising  that  the  potential  at  the  center  of  the 
ring  are  approximately  cylindrical. 

For  low  energy  rings  the  small  toroidal  term  could  be  important  and  could 
have  a  profound  effect  on  the  shape  of  the  orbits.  However,  when  y  »  1,  the 
potentials  for  nQ  —  constant  and  Jg  ■  constant  become  approximately  equal  and 
hence  they  do  not  contribute  substantially  in  Eq.  (7c). 


IV.  Transverse  Ring  Orbits 


Equation  (7b)  has  been  solved  numerically,  using  the  potentials  of  Eq. 
(11).  Typical  macroscopic  beam  orbits  in  the  r,z  plane  are  shown  in  Figs.  6 
to  8.  The  various  parameters  for  those  runs  are  listed  in  Table  II.  Only 
orbits  that  are  at  least  one  beam  minor  radius  away  from  the  wall  are  shown. 
Each  orbit  corresponds  to  a  different  value  of  the  constant  in  Eq.  (7b).  A 
striking  feature  of  the  results  is  the  sensitivity  of  the  orbits  to  the  value 
of  the  constant  in  Eqs.  (7). 

The  number  marked  in  every  fourth  orbit  is  equal  to  10^  [constant  - 
<constant>],  where  the  average  value  of  the  constant,  i.e.  <constant>  for  each 
run  is  shown  at  the  top  of  the  figure.  For  all  the  cases  tested,  less  than  3% 
change  in  the  constant  of  the  motion  was  sufficient  to  generate  orbits  that 
extend  over  the  entire  minor  cross-section  of  the  torus.  Orbits  shown  with 
solid  lines  correspond  to  a  constant  that  is  greater  than  <constant>  and  those 
shown  with  a  dashed  line  correspond  to  a  constant  that  is  less  than 
<constant>. 

All  the  orbits  close  inside  the  vacuum  chamber.  However,  a  fraction  of 
them  lie  inside  the  annular  region  that  extends  from  the  dotted-dashed  line  to 
the  wall.  This  region  has  a  width  that  is  less  than  the  beam  radius  and  hence 
part  of  the  beam  will  strike  the  wall. 

Ring  orbits  in  the  r,z  plane  from  Eq.  (7b)  using  the  potentials  of  Eq. 
(10),  i.e.,  emitting  the  toroidal  terms,  are  shown  in  Figs.  9  to  11.  It  u 
apparent  that  there  is  not  any  noticeable  difference  between  these  orbits  and 
those  of  Figs.  6  to  8. 


The  predictions  of  Eq.  (7b)  are  in  very  good  agreement  with  the  results 
from  our  particle  in  cell  (PIC)  computer  simulation.  Figures  12  to  14  show 
three  computer  simulation  runs.  As  may  be  seen  from  Table  III,  with  the 
exception  of  the  betatron  field,  the  various  parameters  in  the  simulation  are 
the  same  with  those  of  Figs.  6  to  8.  The  slightly  lower  value  of  the  betatron 
field  in  Figs.  12  to  14  is  related  to  the  different  radial  profiles  for  Jq  in 
the  simulation  and  the  potentials  of  Eq.  (11).  The  orbit  wiggles  are  due  to 
the  finite  ring  emlttance,  which  was  taken  zero  in  the  derivation  of  Eq. 
(7b).  It  should  be  noticed  that  in  these  computer  simulation  runs  the 
electron  ring  was  reasonably  well  matched  to  the  magnetic  field  as  it  is 
manifested  from  the  small  variations  in  the  axial  and  radial  ring  envelopes 
shown  in  Fig.  15. 

In  the  general  case,  it  is  difficult  to  derive  an  explicit  expression  for 

the  ring  orbits  in  the  transverse  plane  from  Eqs.  (7b)  and  (11).  However,  in 
2 

the  limit  y  »1,  0Q/  3  3  1  and  v/y  «  1»  such  an  expression  can  be  obtained 
near  the  minor  axis  of  the  torus. 

Since  So  3  3  and  yS  *  y  -  l/2y,  Eqs.  (5)and  (6)  give 


-Uj;  A*Xt  +  -1^  (A®elf  ”  *)  +  •  constant  -  G. 


Expanding  y  near  r  and  using  Eq.  (5),  it  is  obtained 


«r  -  Y-r0  --I4  If  lr  “  +f?-|r 


where  Ar  ■  R-r  .  It  is  shown  in  the  next  section  that  -z — 

o  3r  1  r 


0  and  thus  the 


above  equation  becomes 


jVgVj 


cn.*k>,i 


jki' 


From  Eqs.  (11a)  and  (11c),  the  difference  in  the  self  potentials  can  be 


written  as 


»rlf  -  ♦  -  a*  M  *  *"  ^  A  ^  !  n-V-  (») 

b  a  8a 


Since 


1—  8g  a  1-8  a  l/2y  and  substituting  5y  from  Eq.  (14)  in  the  expansion 


for  l/y  ,  it  is  obtained 


»-»e  *7T  t1 

2\ 


■Af i. 


Similarly,  expanding  l/2y  as 


1  _  1  I  e  I  3*i 

3t  5rT  7TT*  IFlr.  4r> 


o  2y  me 
o 


and  1/R  as 


{■  (■+)  [  i  -f  +  (f)  ] . 


Ar  .  /-Ar - 


o  o 


and  using  a  linear  expression  for  the  external  vector  potential 


o  o 


Eqs.  (13)  to  (19)  give 


p.  tl‘xt  r  vr 

r _ §  go  o  , ,  _ o 

IjJT +  ^  (l  o)  ^ 


vr  n“'r  vr  2  2 

o  i  ( Ar >,  .  r  so  o  o  1  rl z  > 

TtJ  l—)  +  L~27 —  °  ■  T7  J  OH 

Y  a  o  y  a  o 

o  o 


-  [— H_  +  -JL  (Ife-  +  in  — )  1  (— )  ■  G  , 

mCro  2Y  2  2a2  V  ‘  V 


o  2yq  2a 


where 


Ar  -  R-ro,  Az  -  Z,  G  -  G  -  [  1/2yq  +  PQ/mcr 


(v/y  *)  (1/2  ♦  in  -ft)  ♦  Ge  r  /  ],  next  -  lei  Bext/  , 

V  O  '  v  T^J  ZO  0  CJ*  ZO  I  I  20  me* 


and  v  is  the  Budker’s  parameter. 


Equation  (20)  describes  the  ring  orbits  near  the  minor  axis,  when  Y  »  1. 


These  orbits  are  centered  around  the  minor  axis  of  the  torus  when  the 


coefficient  of  the  (—■)  term  is  zero,  i.e.,  when 


r. 

v  u  b 


"cro  2Y  2 

o 


0  9 

For  (rh/_)z  «  1  and  y_  »  1,  Eq.  (21)  predicts  that  — —  «  0.  Therefore,  the 
o  a  o  mcr 


orbits  are  circular  when  the  external  field  index  is  approximately  equal  to  0.5, 


in  agreement  with  the  computer  results  shown  in  Figs.  7  and  13. 


This  result  is  not  in  agreement  with  previous  work,^*^®*^»  which  for  J, 


constant  predicts  circular  orbits  when 


"  "2 


,  1 1  -  <wv  *n  pfr  1 


1  +  (2v/y0)  ( 1/2  +  in  -jS] 

b 


Ic  has  been  determined  that  the  discrepancy  Is  due  to  an  Inconsistency  in 


the  expansion  that  gave  erroneous  results  for  the  two  slow  frequencies 

u>  and  u  and  made  the  expression  for  the  field  index,  i.e.,  Eq.  (22), 

invalid. 

Additional  results  from  Eq.  (7b)  are  shown  in  Figs.  (16)  to  (19).  The 
various  parameters  for  these  runs  are  listed  in  Table  IV.  As  the  ring  current 
increases,  the  orbits  are  dramatically  modified  as  manifested  by  the  results 
of  Fig.  18.  Midway  to  the  wall,  the  orbits  change  from  circles  to  finite 
width  C  shaped  forms  that  evolve  to  cresents  or  ‘’bananas".  At  the  tips  of  the 
crescents  the  bounce  frequency  becomes  zero  and  the  macroscopic  beam  motion 
transitions  from  diamagnetic  to  paramagnetic  and  vice  versa.  Particle  in  cell 
computer  simulation  results  show  that  the  beam  can  go  through  such  a 
transition  without  any  noticeable  interruption.  A  typical  case  is  shown  in 
Fig.  20.  The  various  parameters  of  this  run,  that  lasted  for  more  than  1 
microsecond,  are  listed  in  Table  V.  Figure  20a  shows  the  orbit  of  the  center 
of  the  ring  in  the  r,z  plane.  The  time  interval  beween  two  successive  dots  is 
20  nsec.  According  to  Figs.  20b  and  20c  the  ring  envelope  changes  only 
slightly  during  the  run. 

V.  Extreme  of  the  Constant  of  the  Motion 

The  extreme  of  Eq.  (5)  furnishes  useful  information  on  the  dynamics  of 
the  ring  in  the  r-z  plane.  First,  we  trill  show  that  this  extreme  is  the 
radial  balance  equation  of  motion  for  the  reference  electron. 

Setting  the  partial  derivative  of  Eq.  (5)  with  respect  to  r  equal  to  zero 


13 


"TT»  T77 


9  2  V 

using  the  relation  Y  *  (1  +  S  Y  )  2  and  Eq.  (6),  we  obtain 

_  ».  ext  . .  self 

i. 


where  we  have  assumed  that  3  -  v/c  is  approximately  equal  to  3g  ■  v0/c' 
Substituting  Eq.  (6)  into  Eq.  (24)  and  using  the  equations 


.ext  ..el 

x  «xt  -  4. 1ZL 

Bz  r  +  3r 


(25a) 


self 

.self  A9 
B  ■  — —  — ■  + 

z  r 


(25b) 


E  .  _  i* 

r  3r  » 


(25c) 


it  is  obtained 


Y  m  —  - 


-  el  [E  +  B*,lr>] 

I  L  r  c  z  z  J 


i.e.,  the  radial  balance  equation.  This  equation  gives  the  equilibrium 
position  of  the  ring,  which  is  located  along  the  er  axis.  At  this  position 
the  reference  electron  at  the  centroid  of  the  ring  moves  only  along  the 
toroidal  direction  i.e,  v  ■  v  ■  0. 

1  Z 

When  the  equilibrium  position  is  at  r  -  rQ,  the  toroidal  velocity  of  the 


reference  electron  can  be  determined  from  Eqs.  (6)  and  (21)  and  is 


With  Che  exception  of  Che  very  small  term  on  Che  numerator,  Eq.  (27)  Is  Che 
same  with  the  expression  reported  previously  (20)  for  beams  with  square 
current  density  profile. 

The  external  magnetic  field  required  to  confine  the  ring  at  r  »  rQ 
can  be  readily  found  from  Eq.  (27).  Omitting  the  small  term  in  the  numerator 
of  Eq.  (27), we  obtain 


where  the  single  particle  magnetic  field  is  B*£ 


(28) 


The  magnetic  field  required  to  maintain  the  beam  at  an  equilibrium 
position  that  is  different  than  rQ  can  also  be  determined  from  the  radial 
balance  equation.  Substituting  Ef  and  B*elf  from  Eqs.  (12a)  and  (I2e)  into  Eq. 
(26),  it  is  obtained 


B 


ext 


b“p{  1  t1/2  +  *°-f  + 

b 


R(R-r  ) 

“7 - Z  + 

a  (y&) 


‘b 

— Z - 2 

8a  (yB) 


]  } 


(29) 


Equation  (29)  has  been  derived  under  the  assumption  that  v  is  not  a  function 
of  R. 

The  predictions  of  Eq.  (29)  are  in  excellent  agreement  with  the  result  of 


the  NRL  computer  simulation  code.  Three  examples  are  shown  in  Table  VI 


VI.  Summary 


The  non-linear  beam  dynamics  in  the  plane  transverse  to  the  toroidal 
magnetic  field  is  studied  using  the  two  constants  of  the  motion,  instead  of  the 
linearized  equations  of  motion.  This  approach  allows  the  beam  orbits  to  be 
determined  over  the  entire  minor  cross-section  of  the  torus  and  not  only  near 
the  minor  axis. 

It  was  found  that  the  orbits  are  sensitive  to  the  external  field  index,  to 
the  value  of  the  constant  of  Eqs.  (7)  and  to  the  beam  current.  The  orbits  in 
the  r,z  plane  always  close  inside  the  vacuum  chamber,  although  often  very  near 
the  perfectly  conducting  wall.  As  a  result,  beam  interruption  will  occur 
whenever  the  electron  ring  moves  along  one  of  these  orbits. 

In  addition,  it  has  been  shown  that  the  extreme  of  Eq.  (5)  provides 
information  on  the  external  magnetic  field  required  to  confine  the  ring  at  its 
equilibrium  position  and  the  displacement  of  the  equilibrium  position  when  the 
beam  energy  is  not  matched  to  the  vertical  field. 


Y' 


JvV'iO 


Table  1 


Parameters  for  the  runs  shown  In  Figs.  (2)  to  (5) 


FIs.  2  Fig.  3 


Parameter 


Torus  major  radius  rQ  (cm) 

Torus  minor  radius  a  (cm) 

Electron  ring  minor  radius  r^  (cm) 

Electron  ring  vertical  displacement  Z  (cm) 


R 


Fig.  5 


100 

100 

— 

32 

16 

16 

16 

3 

3 

3 

n 


Fig.  6 


External  field  index 

0.35 

0.5 

0.65 

Torus  major  radius  (m) 

100 

100 

100 

Torus  minor  radius  (cm) 

16 

16 

16 

Ring  minor  radius  (cm) 

3 

3 

3 

Ring  current  (kA) 

5 

5 

5 

Electron  energy  (MeV) 

3.123 

3.123 

3.123 

Betatron  field  B  _  (g) 
zo 

138.4 

138.5 

138.5 

Vacro 

-0.0018 

-0.0020 

-0.0023 

Table  III 


Parameters  of  the  computer  simulation  runs  In  Figs.  (12)  to  (14) 


External  field  Index  n 


Initial  beam  energy  YQ  *  7.117 


Beam  current  I  (kA)  ■  5  KA 


Torus  major  radius  r  (cm)  ■  100 


Initial  beam  minor  radius  r^  (cm)  -  3 


Torus  minor  radius  a  (cm)  ■  16 


Betatron  magn.  field  at  r  ,  z  •  o,  B  (G)  *  136.2 


Toroidal  magn.  field  at  r„,  z  ■  o,  B  .  (G)  ■  388 


Initial  emlttance  e  (rad  -  cm)  >0.1 


Table  IV 


Parameters  for  the  runs  shown  in  Figs.  (16)  to  (19) 
Cylindrical  Potentials 


Parameter 

■HUH 

Pig.  18 

fiBBli 

Hg  It 

Torus  major  radius  rQ  (cm) 

100 

100 

100 

Torus  minor  radius  a  (cm) 

16 

16 

16 

Electron  ring  minor  radius  r^  (cm) 

3 

3 

3 

Electron  ring  energy  E  (MeV) 

1.0 

1.0 

3.0 

Electron  ring  current  I  (kA) 

1.0 

2.0 

10.0 

External  field  index  n 

0.5 

0.5 

0.5 

Equilibrium  position 

100 

100 

100 

Vertical  magnetic  field  B  (G) 

zo 

51.74 

56.08 

159.24 

Norm,  canonical  angular  momentum  Pg^>cr 

0 

0 

0 

0 

Initial  beam  energy  yQ  -  2.76  (E  ■  0.9  MeV) 

Beam  Current  1  (kA)  -  1 

Major  radius  rQ  (ca)  -  100 

Initial  beam  minor  radius  r^  (cm)  >2.5 

Torus  minor  radius  a  (cm)  -  16 

Initial  beam  center  position  r^  (cm)  *  111.0 

Betatron  magn.  field  at  rQ,  z  -  o,  BQZ  (G)  -  47 

Toroidal  magn.  field  at  rQ,  z  -  o,  Bofl  m  400 

Initial  emittance  e  (rad  -  cm)  ■  0.175 

Initial  temperature  spread  (half-width)  — ^  ■  0 

External  field  index  n  *  0.5 

Self  field  index  ng  -  8.6 

Wall  conductivity  -  • 

Time  step  (nsec)  ■  0.10 


No.  of  particles  ■  2048 


Perfectly 

Conducting 

Torus 


-  [A*e'f(R,Z)  /2|e|N  0a] 


^  f,  <r~  *  - 


x  Computer 
-  Analytical 


R  (cm) 


Fig.  2.  Normalized  self  magnetic  vector  potential  from  Eq.  (10a)  and  from 
computer  code  PANDIRA.  The  various  parameters  are  listed  in 


Table  I. 
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Normalized  self  electric  potential  as  a  function  of  R  from  Eq.  (10b)  and 
computer  code  PANDIRA.  The  various  parameters  are  listed  In  Table  I. 
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AVERAGE  BEAM  POSITION 


Computer  simulation  results  showing  the  orbit  of  the  ring  centroid  In  the 
transverse  plane  for  n  «  0.35.  The  various  parameters  for  this  run  are 
listed  In  Table  III.  The  Initial  ring  position  was  R  -  108  cm.  Z  »  0. 


Constant  of  the  Motion  -  <  3.492  > 


Fig. 17.  Orbits  of  the  ring  centroid  as  in  Fig.  16,  but  with  1*2  KA. 


Constant  of  the  Motion  -  <  9.457  > 


Constant  of  the  Motion  -  <  3.103  > 
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C0NT0U9  FROM  -0.48000E-01  T0  0.20000E-01  CONTOUR  INTERVAL  V  0.40000E-02  RT(3.3)»  0.411186-02  LABELS 

SCALEO  BY  10000. 


19.  Orbits  of  the  ring  centroid  as  in  Pig.  16,  but  with - -  *  0.  The 

ncr 

equilibrium  position  is  shifted  to  the  right.  ° 


Fig. 20a  Computer  simulation  results  showing  the  orbit  of  the  ring  centroid  In  the 
transverse  plane.  The  time  Interval  between  two  dots  is  20  nsec.  The 

various  parameters  for  this  run  are  listed  in  Table  V.  About  40  nsec 
before  the  end  of  the  run  the  bounce  frequency  became  zero  without  any 
noticeable  disruption  of  the  ring. 


Vf,' 


iiilM 


Fig. 20b  Radial  ring  envelope  for  the  run  of  Fig.  20a. 
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